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Abstract 

We compute the full set of weak-scale gauge and Yukawa threshold corrections 
in the minimal supersymmetric standard model, including all finite (non-logarithmic) 
corrections, which we show to be important. We use our results to examine the effects 
of unification-scale threshold corrections in the minimal and missing-doublet SU(5) 
models. We work in the context of a unified mass spectrum, with scalar mass Mq 
and gaugino mass M1/2, and find that in minimal SU(5) with squark masses less than 
one TeV, successful gauge and Yukawa coupling unification requires *S ^0 and 

Mq — 1 TeV. In contrast, we find that the missing-doublet model permits gauge and 
Yukawa unification for a wide range of supersymmetric masses. 



1 Introduction 



With the advent of precision measurements at LEP and the observation that the gauge cou- 
phngs unify in the minimal supersymmetric standard model there has been a resurgence 
of interest in supersymmetric grand unified theories. In most analyses the values of the elec- 
tromagnetic coupling, OEM, and the standard-model weak mixing angle, s^j^ = sin^ 9q,m{Mz), 
are taken from experiment, and converted into supersymmetric DR parameters, a and 
= sin^ 9{Mz), using the leading-logarithmic contribution to the supersymmetric threshold 
corrections. The DR parameters are then used to determine the MS strong coupling constant 
at the Z-scale, as{Mz), ignoring all unification-scale threshold corrections. Alternatively, 
the measured value of as{Mz) is used to constrain the unification-scale parameter space. 

In this letter we will take a closer look at this procedure, and report on the results of a 
complete next-to-leading-order analysis of supersymmetric unification. Our approach is new 
in two respects. First, we include all finite corrections. This implies that our weak mixing 
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angle is related to Sgj^ as follows 



s = Sgj^ + leading log + finite, (1) 

where the leading logarithms are of the form log(MsusY/^z) and the finite corrections are 
of order M|/MgugY- If all the superpartner masses are heavier than a few times Mz, the 
finite corrections are negligible, in accord with the decoupling theorem. However, realistic 
supersymmetric models typically have light particles with masses near Mz, so the finite 
corrections can be significant. 

Second, we do not use a combined-fit value of Sg^j to compute because the finite term 
in eq. (|l|) is different for each observable. Therefore, in our analysis, we calculate directly 
from a given set of inputs, namely, the Fermi constant Gp, the Z-boson mass Mz, the 
electromagnetic coupling a em, the top-quark mass m^, and the parameters that describe 
the supersymmetric model. (See also ref. 0.) 

We emphasize that a precise evaluation of is essential because the renormalization 
group equations (RGB's) can naturally amplify small corrections. For example, the one- 
loop renormalization group equations, dcji/dt = /Sj^f/lGvr^, together with gauge coupling 
unification, imply 

g'M giif^) g'M ' ^ ^ 

at any scale /i. Solving for as{Mz) and varying the inputs a and s^, we find 



(5a, 6a 6s 
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— - 7.5 — . (3) 



as{Mz) a s 

This shows that a 1% error in the determination of leads to an error in the evaluation of 
as{Mz) of 7.5%. 

In what follows, we will briefly describe our calculation of a and s^. We will then 
examine our prediction for as{Mz) in the supersymmetric parameter space. We will find the 
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unification-scale threshold corrections that are necessary for gauge coupling unification with 
a given supersymmetric spectrum, and compare our results with the threshold corrections 
that arise in the minimal and missing-doublet SU(5) models. We will close with a similar 
analysis under the additional assumption of bottom-tau Yukawa unification. 



2 Gauge Coupling Unification 



We start by sketching our calculation of the DR gauge parameters gi{Mz) and g2{Mz)- As 
discussed above, we take as inputs the Fermi constant Gp = 1.16639 x 10^^ GeV^^, the Z- 
boson mass Mz = 91.187 GeV, the electromagnetic coupling a^u = 1/137.036, the top-quark 
mass rrit, and the supersymmetric parameters. From these we calculate the electromagnetic 
coupling a and the weak mixing angle in the DR renormalization scheme, 

" " 1-Ad' ~ y2G^M|(l-Af) ' ^ ' 

where = cos^ 0{Mz), 

Aa = 0.0685 ± 0.0007 - ^j-Ilogf^) + '-^log(^) + ^ log (5) 



and 

Ar = Aa + — — 2 — — — r-^ — + vertex + box . (6) 
M-^y Mz 

Equation (^) includes the light quark contribution extracted from experimental data [Q, 
together with the leptonic contribution. It also contains the logarithms of the 14^-boson, 
top-quark, charged-Higgs, squark, slepton, and chargino masses. In eq. (H), the fl denote 
the real and transverse parts of the gauge boson self-energies, evaluated in the DR scheme. 
Equation also includes the vertex and box contributions that renormalize the Fermi 
constant, as well as the leading higher-order and QCD standard-model corrections given 
in ref . ^ . (For a complete description of our calculation see ref . . Note that in this letter, 
all hatted objects are DR quantities, and all masses are pole masses.) 

From these results we find the weak gauge couplings gi{Mz) and g2{Mz) using the 
DR relations 

g,{Mz) = UMz) = ?, (7) 

V 3 c s 

where a = e'^/Air. These couplings serve as the weak-scale boundary conditions for the two- 
loop supersymmetric RGE's. They implicitly determine the unification scale Mqut through 
the condition ^i(Mgut) = ^2(Mgut) = 5'gut- 
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We fix the strong coupling ^3(Mgut) at the scale Mqut by the unification condition 

^3(Mgut) = 5'gut(1 + eg) , (8) 

where eg parametrizes the mo del- dependent unification-scale gauge threshold correction. We 
then run ^3(Mgut) back to the Z-scale to find the MS coupling 

a,(Afe) = (Ml^) (1 - A„.)- , (9) 

where Aa^ denotes the weak-scale threshold correction 

Aa. = ^P + ^logf^Ui:llogf!^U21ogf^)) . (10) 

This procedure allows us to determine as{Mz) for a fixed eg and a given supersymmetric 
spectrum. In what follows, we will make the additional assumption that the supersymmetric 
spectrum unifies as well. Therefore we also assume that the three gaugino masses unify to 
a common value of M1/2 = Mi/2(Mgut), the scalar masses unify to a common scalar mass 
Mq = Mo(Mgut)5 and the soft trilinear scalar coupling parameters unify to Aq = Ao(Mgut)- 

We evolve the parameters M1/2, Mq and Aq to the weak scale using the two-loop super- 
symmetric RGB's 0. We require the parameters to be such that electroweak symmetry is 
spontaneously broken, as naturally occurs when the top-quark mass is large. Then the Higgs 
bosons Hi and H2 acquire expectation values Vi = Vi{Mz) and V2 = V2{Mz); their ratio is 
denoted tan/? = tan/3(Mz) = V2/V1. 

We extract the supersymmetric masses from the running DR parameters and tan/3 at 
the scale Mz- (We choose /i > 0, where the superpotential contains the term +fieijHlH2, 
with ei2 = +1.) We then substitute these masses into eqs. (4) - (6), and repeat the entire 
procedure to find a self-consistent solution to the renormahzation group equations. Typically, 
the process converges after only a few iterations. It allows us to predict as{Mz) consistently, 
including all finite corrections. 

As a point of reference, we show in Fig. 1 our results for as{Mz) in the Mq, M1/2 plane, 
in the absence of unification threshold corrections, for rrit = 170 GeV, tan [3 = 2, and Aq = 0. 
Comparing with the experimental value as{Mz) = 0.117 ± .005 we see that as{Mz) is 
rather large.Q If, for naturalness, we require the squark masses to be below 1 TeV, we obtain 
the lower bound 

asiMz) > 0.126 (No unification thresholds, < 1 TeV, = 170 GeV), 

assuming the validity of perturbation theory, which we take to mean that the DR top-quark 
Yukawa coupling A((Mgut) < 3. If we tighten this condition to A((Mgut) < 1, the bound 
increases by 0.002. For smaller rrit, as{Mz) is smaller. For example, if rrit = 160 GeV, 
the limit reduces by 0.002. These variations apply independently of the unification-scale 
thresholds. Our numbers agree qualitatively with the results of ref. 0]. 

"'^The experimental error in Aa leads to a ±0.001 uncertainty in as{Mz). 
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Figure 1: Contours of as{Mz) in the Mq, M1/2 plane with rrit = 170 GeV, tan P = 2 
and ^0 = 0. The region ruq < 1 TeV is bounded by M1/2 ^ 400 GeV and Mq < 1 
TeV. 



The values of as{Mz) quoted here are larger than in many previous analyses for two 
reasons. First, during the past few years the central value of Sgj^, as determined from 
precision electroweak measurements, has been decreasing. (This change is correlated with 
the increasing best-fit value of rrit.) A smaller value of Sg^j leads to an increase in as{Mz), 
as can be seen from eq. (^. Second, the finite corrections decrease s^, and therefore increase 
as{Mz)- The finite corrections are important in the region M1/2 ^200 GeV where as{Mz) 
is appreciably larger than in the leading logarithmic approximation, as shown in Fig. 2. 

Of course, the value of as{Mz) can be reduced by a unification-scale threshold correction 
with tg < 0. In the minimal SU(5) model 0, the unification-scale gauge threshold correction 
is 

where M^^ is the mass of the color-triplet Higgs fermion that mediates nucleon decay.0 From 
this expression, we see that < whenever Mh^ < Mqut- However, M^^, is bounded from 
below by proton decay experiments, so > throughout most of the parameter spacej^ 

^The gauge threshold e'g also receives contributions from higher-dimensional Planck-scale operators that 
we ignore. They lead to a ±0.006 uncertainty in as{Mz) [0. 

•^We use the formulae given in ref. and choose the conservative values f3 — 0.003 GeV'^ and = 
0.4. 
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Figure 2: The MS coupling as{Mz) vs. Mi/2- The curve labeled LLA shows the 
result if we include only the logarithms of the supersymmetric masses (the leading 
logarithm approximation) , while the solid line corresponds to the full result including 
all finite corrections. 



In minimal SU(5) we find that the smallest possible values for as{Mz) are typically even 
larger than those in Fig. 1. The only exception occurs for Mq ^ M1/2, where the proton 
decay amplitude is suppressed. This determines the lower limit 

asiMz) > 0.123 (Minimal SU(5), < 1 TeV, rrit = 170 GeV), 

assuming Ai(MGUT) < 3. In fact, as long as rrig < 1 TeV, as{Mz) < 0.127 can only be 
obtained in the region Mq ~ 1 TeV. For example, if Mq < 500 GeV, as{Mz) > 0.128. 

The missing-doublet model is an alternative SU(5) theory in which the heavy color- 



triplet Higgs particles are split naturally from the light Higgs doublets [12]. In this model 



the unification-scale gauge threshold can be written as W4 



Thus, for fixed Mfj^, the missing-doublet model has the same threshold correction as the 
minimal SU(5) model, minus 4%. In eq. Mf^ is the effective mass that enters into the 



proton decay amplitude, so the bounds on Mh-^ in the minimal SU(5) model also apply to 
Mf^ in the missing-doublet model. 

The large negative correction in eq. (p!2D gives rise to much smaller values for as{Mz)- 
This is illustrated in Fig. 3, where we show the upper and lower bounds on eg in the min- 
imal and missing-doublet SU(5) models, together with the values of eg necessary to obtain 



5 



0.04 



0.02 



0.00 



Mo=100 GeV 



-0.02 



0.04 



iiiiiiiiiiiiiiiiiiiiiiiii 



I I 1 1 1 1 1 1 1 1 1 1 1 1 1 iiiiiiiiiiiiiiiiiiiiiii 



(a) 



minimal SU(5) 




missing doublet SU(5) 




50 100 200 500 

Mwg (GeV; 



Mw2=100 GeV 



(b) 



minimal SU(5) 




missing doublet SU(5) 




50 100 200 500 1000 

Mo (GeV) 



Figure 3: The light shaded regions indicate the allowed values of the gauge coupling 
threshold correction Cg in the minimal and missing-doublet SU(5) models. The dark 
shaded region indicates the range of €g necessary to obtain as{Mz) = 0.117 it 0.1. 



as{Mz) = 0.117 ± 0.01. For both SU(5) models we bound eg from below by the limits on 
proton decay and from above by the requirement Mh^ < 10^^ GeV. Note that for the missing- 
doublet model, the region of allowed eg nearly overlaps the region with as{Mz) = 0.117±0.01, 
regardless of the supersymmetric particle masses. 

3 Yukawa Coupling Unification 

In typical SU(5) models the bottom and tau Yukawa couplings unify at the scale Mgut- 
In the remainder of this letter, we will investigate bottom-tau unification to the same next- 
to-leading-order accuracy that we used in our analysis of the gauge couplings. Our anal- 
ysis is based on two-loop RGB's and on the complete one-loop expressions that relate the 
DR Yukawa couplings to the physical fermion masses 0. 

Let us illustrate our procedure for the case of the tau Yukawa coupling. We start by 
determining the DR expectation value v from the Z-boson mass, 

^(g'^ + f)v^ = M| + nz(M|) . (13) 

Then, given the pole mass mr = 1.777 GeV [|], we find the tau Yukawa coupling A^(M^) 
using the DR relation 

Xr{Mz) V cos /3/V2 = mr + trinir) , (14) 
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where ST-(mT-) = Si +mT-S^, and the tau self-energy is Ei +|^S^ + 75 (S5 + We follow 

a similar procedure to find the top-quark Yukawa coupling Ai(M^). 

Once we have the top and tau Yukawa couplings at the scale M^, we evolve them to the 
scale Mgut using the two-loop supersymmetric RGE's 0. At Mgut we fix the DR bottom- 
quark Yukawa coupling ^^(Mgut) through the unification condition 

Af,(MGUT) = A,(MGUT)(l + eb) , (15) 

where parametrizes the unification-scale Yukawa threshold correction. We then run all 
the Yukawa couplings back to the Z scale, and iterate the procedure self-consistently to 
determine Xb{Mz)- Finally, we apply the weak-scale threshold corrections to find the pole 
mass for the bottom quark. 

For a top-quark pole mass rrit = 170 GeV, this procedure typically gives a bottom-quark 
mass outside the range of experiment (which we take to be 4.7 < nib < 5.2 GeV [§). From 
previous analyses we know that for nit < 200 GeV there are two regions of tan/3 where 
bottom-tau unification might occur: tan/3 < 2 and tan/3 > 40. Here we focus on the 
small tan /? branch, where the top-quark Yukawa coupling at the unification scale is large. 
At (Mgut) > 1- 

As a point of reference, we first present our results with no unification-scale threshold 
corrections. In Fig. 4 we show nib and as{Mz) versus mt, for various values of tan/3, Mq, 
and Ml/2, with A((Mgut) = 3. From the figure we see that the bottom-quark pole mass is 
generally too large, unless the squark masses are of order 1 TeV. 

The value ^^(Mgut) = 3 was chosen because it lies on the edge of perturbation theory. 
Smaller values of Ai(MGUT) give rise to larger values of nib, so the curves in Fig. 4 can be 
interpreted as lower limits on the bottom-quark mass. If we require the squark masses to be 
below 1 TeV, we obtain the lower bound 

nib > 5.1 GeV (No unification thresholds, nig < 1 TeV, nit = 170 GeV). 

In the absence of unification- scale threshold corrections, we see that, in the small tan/3 
region, bottom-tau Yukawa unification favors large values of and large values of the 
supersymmetric mass scale. The above bound depends slightly on the top-quark mass, and 
sensitively on At(MGUT)- For mt = 160 GeV, the limit reduces by 0.1 GeV. The bound 
increases by 0.3 GeV for ^^(Mgut) = 2 and by 0.8 GeV for Af(MGUT) = 1- These variations 
also apply for 7^ 0. 

As with as{Mz), the picture is altered by unification-scale threshold corrections. To 
understand their effects, we first note the striking similarity between the nib and as{Mz) 
curves in Fig. 4, which implies that the value of rrib is tightly correlated with the value of 
as{Mz). This leads us to expect that the gauge threshold correction will have an important 
effect on nib- 

In Fig. 5 we show the prediction for versus mt, for fixed = 0, with Af(MGUT) = 2, 3 
and eg chosen to give a fixed value of as{Mz). We see that any model which gives as{Mz) — 
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Figure 4: The bottom-quark mass and as{Mz) vs. mt for the case of no unification- 
scale thresholds, for various values of tan/?, with Aq = and A4(Mgut) = 3. The 
right (solid) leg in each pair of lines corresponds to M1/2 varying from 60 to 1000 
GeV, with Mq fixed at 60 GeV. The left (dashed) leg corresponds to Mq varying from 
60 to 1000 GeV, with Mi^2 — 100 GeV. On the solid lines the circles mark, from 
top to bottom, Ml/2 = 60, 100, 200, 400, and 1000 GeV, and on the dashed hues 
the circles mark Mq = 60, 200, 400, and 1000 GeV. Note that the lowest point on 
each left leg and the second-to-lowest point on each right leg corresponds to rriq ~ 1 
TeV. The horizontal dashed lines indicate rrib = 5.2 GeV and as{Mz) = 0.127. The 
X 's mark points with one-loop Higgs mass ruh < 60 GeV. 



0.12 also predicts mf, ~ 5 GeV, with = 0. Furthermore, any —23% < e?, < 9% gives an 
acceptable value for nih in this case. Alternatively, if as{Mz) — 0.13, a large and negative 
is required to achieve the central value for nib, namely —33% < e;, < —8%. Figure 5 shows 
that these conclusions hold independently of the top-quark mass and the supersymmetric 
mass scale. 

In Fig. 6 we plot the most favorable values of as{Mz) and nib in the minimal SU(5) 
model. In this model the Yukawa threshold is given by |15| 

where My is the mass of the superheavy SU(5) gauge bosons. We define the most favorable 
value as follows. We first minimize the value of as{Mz) by picking the smallest possible 
Mh-j consistent with nucleon decay. We then minimize nib by choosing the smallest My 
allowed by the constraint that the SU(5) model Yukawa couplings remain perturbative, 
My > max(0.3M//3, O.SMgut)- From the figure we see that this brings nib to an acceptable 



8 



> 



6.0 



5.5 



5.0 



4.5 



4.0 



tailjS = 1.3 1.5 2.0 2.5 

/a=0.13" 

e e e g \ ^ 

Q e B e 7 

/a =0.12 

{a=o.ir 



160 



180 



200 



220 



GeV 



Figure 5: The bottom-quark mass vs. the top-quark mass for fixed values of as{Mz) 
and various tan/3. The solid lines correspond to At(MGUT) = 3 while the dashed 
lines correspond to Af(MGUT) = 2. The upper line in each pair corresponds to a 
light supersymmetric spectrum with Mq = M1/2 = 80 GeV. The lower line in each 
pair corresponds to a heavy spectrum, Mq = 1000 GeV, M1/2 = 500 GeV. 



range for squark masses of order 1 TeV. Indeed, we find the limitQ 

nib > 5.1 GeV (Minimal SU(5), nig < 1 TeV, rrit = 170 GeV). (17) 

In contrast to minimal SU(5), the missing-doublet model can readily accommodate 
bottom-quark masses in the range 4.7 GeV < rrih < 5.2 GeV, with e^, ~ 0. For exam- 
ple, with rrit = 170 GeV and Mq = M1/2 = 100 GeV, minimal SU(5) requires ef, in the 
range —15 to —65%. The missing-doublet model, however, requires e?, in the range —30 to 
+20%. Hence even for small (9(100 GeV) supersymmetric masses the missing-doublet model 
naturally accommodates both gauge and Yukawa coupling unification. 

4 Conclusion 

In this paper we have computed the complete one-loop weak-scale threshold corrections in 
the minimal supersymmetric standard model. We used them to study gauge and Yukawa 
unification with and without unification-scale threshold corrections. In the absence of such 
corrections we find that as{Mz) and rrif, are large unless the squark masses are larger than 

"^We could try to decrease further by increasing Mh^, which decreases e'j,. However, the corresponding 
increase in as{Mz) compensates for this so that the limit in eq. remains unchanged. 
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Figure 6: The same as Fig. 4, with the most favorable minimal SU(5) threshold 
corrections, as defined in the text. 



about 1 TeV. Adding minimal SU(5) threshold corrections, and requiring < 1 TeV, we 
find that a,{Mz) < 0.127 occurs only for M1/2 < Mq and Mq ~ 1 TeV. Additionally, the 
condition nif, < 5.2 GeV is fulfilled only for rrig > 1 TeV. In the missing-doublet model, 
however, the threshold corrections permit acceptable values for as{Mz) and rrib for either a 
light or heavy supersymmetric spectrum. 
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